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We investigate the influence of the electron-phonon interaction on the decay dynamics of a quan- 
tum dot coupled to an optical microcavity. We show that the electron-phonon interaction has 
important consequences on the dynamics, especially when the quantum dot and cavity are tuned 
out of resonance, in which case the phonons may add or remove energy leading to an effective non- 
resonant coupling between quantum dot and cavity. The system is investigated using two different 
theoretical approaches: (i) a second-order expansion in the bare phonon coupling constant, and 
(ii) an expansion in a polaron-photon coupling constant, arising from the polaron transformation 
which allows an accurate description at high temperatures. In the low temperature regime we find 
excellent agreement between the two approaches. An extensive study of the quantum dot decay 
dynamics is performed, where important parameter dependencies are covered. We find that in 
general the electron-phonon interaction gives rise to a greatly increased bandwidth of the coupling 
between quantum dot and cavity. At low temperature an asymmetry in the quantum dot decay 
rate is observed, leading to a faster decay when the quantum dot has a larger energy than to the 
cavity. We explain this as due to the absence of phonon absorption processes. Furthermore, we de- 
rive approximate analytical expressions for the quantum dot decay rate, applicable when the cavity 
can be adiabatically eliminated. The expressions lead to a clear interpretation of the physics and 
emphasizes the important role played by the effective phonon density, describing the availability of 
phonons for scattering, in quantum dot decay dynamics. Based on the analytical expressions we 
present the parameter regimes where phonon effects are expected to be important. Also, we include 
all technical developments in appendices. 

PACS numbers: 78.67.Hc, 03.65.Yz, 42.50.Pq 



I. INTRODUCTION 

The study of cavity QED (cQED) has for decades been 
an important topic in physics. Originally, the main in- 
gredients were atoms, highly confined modes of light, 
and their mutual interaction. Recent years have seen 
a rebirth of cQED, but with focus shifted from the pure 
setting of atoms and cavities, to the complex setting of 
many-body physics found in semiconductor solid-state 
systems. A major driving force behind this shift is the 
advent of quantum information technologies^, with the 
requirements of applications pushing for the exploration 
of new material platforms. A scalable all-solid-state plat- 
form, where the interaction between light and matter can 
be engineered and controlled to a high degree^^, could 
help usher practical devices employing quantum informa- 
tion technologies. 

A solid-state platform, however, also poses new chal- 
lenges owing to its inherent many-body nature, namely 
the effect of the environment on the fragile quantum 
states of light and matter and their coherent interaction, 
which are essential for many applications. Several recent 
studies®^ have shown that simple concepts useful in un- 
derstanding atomic cQED systems break down on both a 
quantitative and qualitative level for all-solid-state cQED 
systems. The two main reasons for the departure from 



the usual picture are (i) the impossibility of quantum 
emitters in the solid-state to be described as simple two- 
level systems and (ii) the stronger coupling to structured 
environments in the form of, e.g., phonons and electronic 
inter-particle Coulomb interactions. 

For an all-solid-state cQED system consisting of a 
semiconductor quantum dot (QD) and an optical micro- 
cavity especially the interaction with phonons has at- 
tracted a considerable amount of attentio n. I t has been 
shown to influence cQED emission spectraP^HUS t gj V e 
rise to detuni ng dependent spectral asymmetries in QD 
lifetime^ 6 * 7 * 11 !, as well as yielding unexpected broadening 
mechanisms in connection with Mollow triplets for coher- 
ently driven systems^HUl. 

The majority of studies has focused on the effect of 
phonons in the spectral domain, where typically the spec- 
trum of the emitted light from the entire cQED system is 
collected and analyzed. However, for quantitative stud- 
ies, measurements in the temporal domain are in many 
cases expected to be superior-^ due to their insensitivity 
towards collection efficiencies. The study of spontaneous 
emission decay has been employed to probe the environ- 
ment in which the emitter is emerged into, be it, e.g., 
electromagnetic^ or plasmonicPSI in nature. 

In a previous study^ we showed how, at low temper- 
atures, the phonon interaction gives to a significantly 



faster decay of an excited QD, whose transition frequency 
is blue-shifted relative to the cavity, as compared to a red- 
shifted QD. In addition, coupling to phonons gives rise to 
a renormalization of the light-matter coupling strength. 
Similar results have independently been obtained by 
other a 6 l 19 l It was argued that the non-trivial phonon ef- 
fects could only be explained if the phonons were treated 
as interacting with the electron-photon quasi-particle, 
the polariton, and not with the bare electron^. 

Here, we present the details of the theory developed 
in Ref. [7] and expand the treatment by comparing to an 
alternative method, more appropriate for higher temper- 
atures. Excellent agreement between the two methods is 
found in the low-temperature regime, which is of our pri- 
mary concern. We perform an extensive parameter study, 
providing, a good picture of the dynamics in different 
regimes. We furthermore derive an analytical expression 
for the QD decay rate, which makes the involved phys- 
ical processes apparent. The analytical expression has 
very recently been used to experimentally map out the 
effective phonon density^. Furthermore, it inspired to a 
novel approach for decreasing phonon-induced dephasing 
in cQED systems^!. Furthermore, we provide a simple 
explanation as to why phonon-induced asymmetries have 
largely remained unobserved in experimental data until 
recently. 

The paper is organized as follows. In Section [TT] we 
describe the model, emphasizing the interaction with 
phonons, and introduce the the polaron transformation 
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enabling the treatment of higher temperatures. Sec 
gives a detailed description of the theoretical formalisms 
employed as well as providing a physical interpretation 
of the resulting equations of motion. We pursue two 
methods; The first is based on a second-order expan- 
sion in the phonon coupling, yielding simple equations 
that provide valuable insight into the physics. The sec- 
ond method employs a partially infinite order expansion 
in the phonon coupling, based on the polaron transfor- 
mation, which leads to more accurate results but less 
physical insight. In Section [IV| we present a detailed pa- 
rameter investigation of the models, covering experimen- 
tally relevant parameter regimes and discuss the physics 
of the system. Furthermore, we perform a large detuning 
expansion and obtain analytical expressions for the total 
decay rate of the QD, which explicitly accounts for the 
different contributions to the system decay and make the 
physical processes very apparent. Finally in Section [V] 
we summarize and conclude. 



II. MODEL SYSTEM 

In this section we present the model used to de- 
scribe the cQED system, including the interaction with 
phonons. The system is illustrated schematically in 
Fig. [T] We also devote a section to the polaron trans- 
formation. 
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FIG. 1. Schematic of the cavity QED system including the 
phonon interaction. The QD-cavity coupling strength is g 
and the QD-phonon interaction matrix elements are M*.. The 
rates V and k yield decay of the QD and cavity, respectively. 
Pure dephasing of the QD is included through 7 and A is the 
QD-cavity detuning. 



A. Cavity QED system 

The part of the system consisting of the QD and cavity 
can be represented by the Hamiltonian 



H, 



cQED 



(1) 



Here H s describes the QD-cavity system and H 1 , H K , and 
Hr describe various interactio ns wit h the environment, 
included as Lindblad loss term^21122l ! to be discussed be- 
low. The QD-cavity Hamiltonian reads 

H s = ^2 hjjic\ci+huj cav a) a + hg(a l c^c e +clc s a), (2) 



where the usual rotating wave and dipole approximations 
have been applied. The energy of the ground (excited) 
QD state is hui g (Huj c ) with corresponding fermionic op- 
erators c|,c g (cj,c c ), the energy of the cavity photon is 
huj cav with corresponding bosonic operators a\a, and 
g is the interaction strength between the cavity pho- 
ton and the electron in the QD. As we are only con- 
cerned with the dynamics of the system on the single 
photon level, it is advantageous to project the second 
quantized Hamiltonian, presented above, onto a lower 
dimensional Hilbert space. An appropriate basis to 
span this part of the total Hilbert space is the follow- 
ing: {|1> = |e, n = 0) , |2) = |g, n = 1) , |3) = |g, n = 0}}, 
where n refers to the number of cavity photons. If we 
project onto this basis and shift to a rotating frame, we 
can write the QD-cavity system Hamiltonian as 



H s = hAan + hg(a 12 + a 2 i), 



(3) 



where A = uj c — uj g — cj cav = uj cg — w cav is the QD-cavity 
detuning and a pq = \p) (q\ is the standard projection op- 
erator. The detailed steps are given in Appendix [B| 

The remaining terms in -Hcqed all give rise to differ- 
ent forms of losses, which we include through the Lind- 
blad formalism often employed in density matrix the- 
ory. The Hamiltonian iJ 7 represents pure dephasing pro- 
cesses, with rate 7, for all transitions connected to the 
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QD, whereas the Hamiltonians H K and Hr account for 
population decay from the cavity and QD by rates k and 
r, respectively^. These rates are taken as parameters 
with experimentally relevant values. 



B. Phonons 

The Hamiltonians involving phonons are given by 
-ffo.ph = ^2 hu k blb k , (4) 

k 

ffc-ph = ( M gg4 c s + M cAc c ) (bl k + b k ), (5) 

k 

where -ffo.ph describes the free phonons and -ff e -ph de- 
scribes the electron-phonon interaction. It sh ould be no- 
ticed that we assume bulk phonon modes^MUl ^he la 
phonon dispersion relation is assumed to be linear in the 
relevant energy range, w k = c s \k\, with c s the speed of 
sound. b k ,b k are the bosonic operators for the phonons. 
The matri x element M„ v in the electron-phonon interac- 
tion iPMHU 

where d the is mass density, c s is the speed of sound in the 
material, V is the phonon quantization volume, D u is the 
deformation potential, and </v( r ) is the electronic wave- 
function for the state involved in the phonon process. We 
neglect the polar coupling to longitudinal optical (LO) 
phonons due to their large energies, ~ 37 meV, com- 
pared to the energies involved in this model, and hence 
very non-resonant nature. Also, we neglect the piezoelec- 
tric coupling to LA phonons, which has been shown to 
have a small effect for the present system^. 

To model the QD wavefunctions, we consider harmonic 
confinement in the direction perpendicular to the growth 
direction^! and infinite potentials in the growth direc- 
tion. This implies wavefunctions for both the ground 
and excited state of the form 

Mr) = 1/2 exp[-(x 2 +^)/(2/^)] 

X C0S(7TZ/Z cff , z ), \z\ < l cS ,z/2 (7) 

where the confinement lengths l xy ,v and l c g tZ can be 
chosen to model a specific system. We choose QD 
and phonon parameters suitable for typical InGaAs 
systems^. 

If we take advantage of the fact that we only consider 
a single electron, i.e. c|c g + c|c c = 1, and project onto 
the basis introduced above, we obtain 

ffc-ph - ^11 Mk ( b -k + h k) = <ruB, (8) 

k 



where we introduced the effective matrix element 

M k = M k c - M g fc g . (9) 
The details are presented in Appendix |B) 

C. The polaron transformation 

We start from the following HamiltoniarP3 

H = hAan + hg(a 12 + a 2 \) 

+ a 11 Y,M k {bl k + b k )+Y,f^kblb k , (10) 

k k 

obtained by combining Eqs. ([3|, Jif, and Q. We then 
apply the polaron transformatiori^ 1 ^ | al " ^, where an op- 
erator O transforms as 

6 = e s Oe- s (11) 

where 

nrk 

S = a n C, C = J2*h(bl k -b k ), A fc = — . (12) 

The idea behind the transformation is to remove the term 
linear in the phonon operators in order to arrive at a 
set of equations that is easier to treat. Physically, the 
transformation shifts the phonon modes according to the 
presence of the electron, determined by the operator an. 
From the exponential nature of the transformation oper- 
ator e s , phonon processes are included to infinite order. 
This has the consequence that multi-phonon effects are 
easily included in the theory, allowing for the description 
of experiments performed at high temperatures. We use 
the bar to signify the transformed frame. The Hamilto- 
nian in the polaron frame becomes 

H = H s i + -f/g'-ph' + H ^h, (13) 

with 

H s , = HAa n + % (A) (a 12 + a 2l ), (14a) 

# S '- P h' = hg{a 12 5X + + a 21 SX_), (14b) 

H , P h = ^2 huj k b k b k . (14c) 

k 

It should be noted that a constant energy shift, induced 
by the phonons, has been absorbed in the QD-cavity de- 
tuning A, see Eq. ( |C8[ ). Also, new phonon related oper- 
ators have been introduced 

SX ± =X ± - (X) , (15) 

X± = e ±c , (16) 

where it holds that 

(A) = (X ± ) . (17) 
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The brackets denote the expectation value with respect 
to the thermal density matrix for the phonons, more pre- 
cisely (• • •) = Tr ph {p p h,o ' ' '}• The detailed derivation 
can be found in Appendix [C] and various relevant prop- 
erties of the operators X± are described in Appendix |E| 
Due to the polaron transformation, the division of the 
total Hamiltonian into a QD-cavity system part and a 
phonon part is no longer possible. Indeed, the new sys- 
tem Hamiltonian, H s i , contains the phonon quantity (X) 
which is seen to renormalize the light-matter coupling 
strength g. It should also be noted that defining the 
new system Hamiltonian in this way, we include pho- 
ton processes to infinite order and respect the detailed 
bala nce c ondition^. From the expression for (X), see 
Eq. (|E14[), it is obvious that < (X) < 1. The in- 



teraction with phonons will thus always decrease the ef- 
fective light-matter coupling strength as a consequence 
of this. The new interaction Hamiltonian, i? s '_ p h', con- 
tains the phonon fluctuation operators 8X±, describing 
fluctuations of the phonon bath around its equilibrium 
value, as well as the light-matter coupling strength g. 



III. EQUATIONS OF MOTION 

In this section we present the theoretical formalism em- 
ployed for analyzing the system, described by the Hamil- 
tonians of the previous section. The explicit form of the 
equations of motion is also presented. 



A. Time-convolutionless approach 

Our basic approach is to set up an equation of mo- 
tion for the reduced density matrix (RDM) of the QD- 
cavity system, where the phonon degree s of free dom are 
traced out. This is a standard techniqu e 1 21 * 22 * 35 ! in which 
the effect of the reservoir enters through various scat- 
tering terms in the equation of motion (EOM) for the 
RDM. These scattering terms can be derived by two dif- 
ferent approaches. In the first, known as the Nakajima- 
Zwanzig projection operator technique^!, the resulting 
EOMs have memory: the present state of the system 
thus depends on the past history. In the second, known 
as the time-convolutionless approach^ (TCL), the EOMs 
are time-local and therefore do not have memory, how- 
ever, the scattering rates become time-dependent. Both 
of these approaches yield, without further approxima- 
tions, a non-Markovian description of the dynamics. 

In this paper we employ the TCL up to second order 
in the perturbation, for the following two reasons: The 
first and most important is that in the limit where the 
light-matter coupling tends to zero, our model reduces 
to the so-called independent boson modeP^. This model 
is known to be exactly solvable using a number of meth- 
ods, one being the second order TCL?' 1 . Even though 
the present model can not be solved exactly using the 
second order TCL, we expect the result to be more ac- 



curate compared to that obtained using the method in- 
volving memory integrals, since that method does not 
lead to the exact solution to second order for g — > 0. 
Other studies have also shown the TCL to be superior to 
the corresponding equation with memory^. The second 
reason is purely practical, in that time-local equations 
are simpler to solve than equations containing memory 
integrals 

The EOMs arising from the TCL may be derived in 
a completely general frameworlP^l, however we follow a 
less rigorous approach in deriving the TCL and present 
the resulting formulas in Appendix [Xj 



B. Phenomenological losses 



As mentioned in Section II A we also include inter 



actions with other reservoirs than phonons to simulate 
a real system with losses. These are included using the 
Lindblad formalism^!] where terms of the form 



L{0, 7 }p(t) = 



- [o ] o P {t) 



p{t)0^0 -20p(t)0^] , (18) 



are added to the EOM for p(t), where p(t) = Tr ph {x(t)} 
is the RDM for the QD-cavity system, x(i) is the den- 
sity matrix for the total system, and Tr p h {• • •} denotes 
the trace operation with respect to the phonon degrees 
of freedom. The above leads to decay with rate 7 of 
the transition corresponding to the operator O. This ex- 
pression may be obtained by taking the white noise, or 
equivalently zero memory, and zero temperature limit of 
the scattering terms presented in Appendix |A"| 

The decay of the cavity field through leaky modes is 
modeled by including the Lindblad term L {032, n} p(t), 
the decay of the excited QD through radiative 
and non-radiative processes is modeled by including 
L {031, r} p(t), and finally a Markovian pure dephasing 
rate is also included through L {en, 27} p(t). We refer to 
Section |II A| for notational remarks. Since LA phonons 
have been included explicitly, and already give rise to 
a pure dephasing rate, it might seem redundant to in- 
troduce an additional pure dephasing channel. However, 
previous work has demonstrated that excited states for 
electrons and holes contribute to pure dephasing pro- 
cesses n ear th e ground state transition energy, due to 
both LA 1221221 and LCf^ phonon interactions. Also, in- 
cluding a finite lifetime of either LO and LA phonons, 
arising, e.g., from anharmonic effect induces a contri- 
bution to the pure dephasing rate^l. For simplicity, we 
assume 7 to be an independent parameter. 



C. Notational remarks 

The resulting EOMs we arrive at are all linear in the 
elements of the RDM. This fact makes it advantageous 
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to formulate the EOMs in the language of linear algebra. 
This can be achieved by mapping the RDM onto a vector 
form as follows 



r W) = [Mt)) , (tr 22 (t)> , (a 12 (t)) , (a 21 (t)) , 

(a 23 (t)),(a 32 (t)),(a 13 (t)),(a 31 (t))] T - (19) 



Here, (<J qp (t)) = Tr s {p(t)<r qp } = p pg (t), where Tr s {---} 
denotes the trace with respect to the QD-cavity basis. 
The QD ground state population, i.e. ((733(4)), has been 
omitted as it does not matter for the dynamics considered 
and may be trivially obtained using the conservation of 
population. The matrix describing the coupling between 
different elements can be divided into three main contri- 
butions 



dt (<r(t)) = [Af coh + M LindWad + M LA (t)] (<x(<)) 

= M(t) (cr(t)) , (20) 



where M co h describes terms originating from the coher- 
ent unitary evolution provided by the QD-cavity Hamil- 
tonian, M Lindblad describes terms from the Lindblad op- 
erators, and Mla(£) describes the time-dependent scat- 
tering terms induced by the coupling to LA phonons. 

As will be shown, M(t) can be written as two decou- 
pled sub-matrices 



M(t) = 



m {ll \t) 
m< 22 )(<) 



(21) 



where m( u )(i) couples the first four elements of (<r(t)), 
to*- 22 -* (t) couples the last four, and all other elements are 
zero. 

In the following two sections we will derive the EOMs 
for the system using the TCL. We present the equations 
arising from the Hamiltonian without the polaron trans- 
formation, denoted the original frame, and with the po- 
laron transformation, denoted the polaron frame. Em- 
ploying the polaron transformed Hamiltonian is expected 
to yield more accurate results compared to the original 
Hamiltonian, especially for elevated temperatures. How- 
ever, the equations resulting from the polaron transfor- 
mation are also more complicated and due to the change 
of basis harder to interpret physically. On the other 
hand, the equations arising in the original frame are sim- 
ple and can be used to gain insight into the physics. 



D. Original frame 

In the original frame, i.e. not employing, the polaron 
transformation, the total Hamiltonian without the Lind- 
blad contributions is 



pii • 



(22) 



where the individual contributions can be found in Eqs. 
([3]), Q, and |8]), respectively. We consider H e _ ph as the 
interaction Hamiltonian, for which the perturbation ex- 
pansion is performed. With this choice only the electron- 
phonon interaction is treated approximately, which is ex- 
pected to be a good approximation, whereas the electron- 
photon interaction is treated exactly and the theory is 
not limited to small values of the light-matter coupling 
strength g. 

To writ e up the TCL EOM for the RDM we use 
Eq. ( Al5[ ) and the time-local scattering term given in 
E q. (A2 0| and_ finally add the Lindblad terms discussed 



in IIIB]to geP 



d tP (t) = -ift- 1 [fl,,p(t)]+5 LA (t) 

+ (L {a 32 , k} + L {a u ,T} + L {a u , 2 7 }) p(t). (23) 

Written in terms of the operator expectation val- 
ues (cr nm (i)), the populations in the QD-cavity sys- 
tem are obtained as follows: The cavity population is 
(a'(t)a(t)) = (<J 22 (t)) and the excited QD population 
(cl(t)c c (t)) = (erii(i)). The off-diagonal elements cor- 
respond to different polarizations or coherences in the 
QD-cavity system, with the relevant one for one-time dy- 
namics being the so-called photon-assisted polarization 
(ci2(i))- Remapping the RDM to vector form, we get 
the following coupling matrices. The coherent terms are 



—ig ig 

ig —ig 

-ig ig iA 

ig —ig — iA. 



(24) 



and 



,(22) 






IfJ 









-ig 



if) 


iA 




' 

-ig 


-iA 



(25) 



and the Lindblad contributions take the form 
diag{M Lindb i ad } = 

- - [2r, 2k, r + k + 2 7 , r + « + 2 7 , 

K,K,L + 2 7 ,r + 2 7 ] , (26) 
and finally the phonon induced terms are 
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- 











-ig>(t) ig<{t) 


-[712(2)- 


U[g>(t)r -i[g<{t)\* 






-[7* 2 (*) + A>oi]J 



(27) 



,(22; 



(*) = 






iQ<(t) 








-i[g<(t)y 






-[713(2) - «A pol ] 










(28) 



polj 



The elements of tola will be defined below, Eqs. (29 1, 
( 30 1 , and ( 32 ) , but first we provide a brief discussion 



of the elements. If one disregards the phonon induced 
scattering terms, these equations constitute the standard 
lossy Jaynes-Cummings model including pure dephasing, 
which has been studied intensely in recent years^HUl. L e t 
us start by discussing the terms in mj^ (t) in more de- 
tail, i.e., the quantities 712(2) and g < (t). If we compare 
the structure of the phonon scattering term, Eq. (27), 
with the non-phonon related terms in the coherent and 
Lindblad contributions to M, a physical interpretation 
of the effects of phonons becomes apparent. 



The rate 712(2) multiplies the photon-assisted polariza- 
tion and therefore the real part of 712 (2) represents pure 
dephasing of this specific polarization, whereas the imag- 
inary part corresponds to an energy shift. The long-time 
limit of this energy shift has been subtracted in the form 
of the quantity A po i = Im {712(00)}, usually referred to 
as the polaron shift, to provide a consistent expansion in 
the electron-phonon interactional. This adjustment has 
been performed everywhere the detuning, A, enters and 
results in an effective QD-cavity detuning close to zero, 
A «Q. 



The quantities g < (t) multiply the populations of the 
excited QD-cavity system in such a way that the real 
part of g < (t) renormalizes the bare light-matter coupling 
strength g. However, in general Re [£/>(£)] 7^ Re [£*<(£)] 
and hence the renormalization does not correspond to an 
overall change in the value of g in the EOM for (1712(2)). 
The imaginary part of C? ^ (i) gives rise to an additional 
decay or growth of the polarization, depending on the 
sign of Im [{/<(£)] , if state 1 or 2 is populated. The influ- 
ence of the degree of excitation in the QD-cavity system 
makes this dephasing channel of a different nature than 
the pure dephasing normally induced by phonons, which 
is well understood, see e.g. Ref. |2"51 



From the scattering term Eq. ( A20 1 we get 

g^(t) = ih- 2 f dt'U^Unit'jD^it'), (29) 
Jo 

712(2) = h- 2 [ dt'[\U 11 (t')\ 2 D<(t') - \U 2l {t')\ 2 D>(t')] 
Jo 



tT 2 f dt'[D<(t') - |C/ 2 i(2')| 2 2Rc[ J D < (<')]] I 
Jo 



A pol = Im {712(00)} , 

7i 3 (2) = h- 2 [ dt'\U n (t')\ 2 D<(t') 
Jo 



h~ 2 f dt'[D<(t') 
Jo 



\U 21 (t')\ 2 D<(t')] 



(30) 
(31) 



(32) 



where it has been used that both D<(t — t') and U nm (t — 
v) only depend on the difference between the two time 
arguments and further the initial time has been assumed 
to be zero. The phonon bath correlation functions enter- 
ing above are defined as 

D*(t) = Y) Mh \ 2 [ n ke ±tuJkt + K + 1) e T ' iWfct ] (33) 
fe 

= ^|M fc | 2 [(2n fc + l)cos(w fc t) Tism{cu k t)} , 
fe 

(34) 

which are related to the phonon bath operators B in the 
following way 



D*(t-1?) = (B(±[t-t , })B(0)) 



(35) 



and rife is the thermal occupation factor for the fc'th 
phonon mode, defined in Eq. (E8 ). The matrix U(t) is the 
time-evolution operator for the QD-cavity system which, 
due to the time-independence of H s , see Eq. can be 
given as a closed form expression 



U(t) = exp(-iH s t/H). 



(36) 



The products of the elements of U(t) occurring in Eqs. 
( 29 1 and ( 30 1 can be interpreted as propagators of the 
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QD-cavity system governed by H Sl representing the pure 
lossless Jaynes-Cummings model. This is easily realized 
by writing the time-evolution of the density matrix for 
the pure Jaynes-Cummings model as 



P 3G (t) = U(t)p' c (0)U\t). 



(37) 



If we assume that p (0) = a kk , i-e. the time-evolution 
starts with the excitation in a single state, we get 

= u nk (t)uljt) = u nk {t)u* mk {t). 



rnm \ 



(38) 



The time-evolution of p^{t,(j kk ) contains the light- 
matter coupling, and so do the Jaynes-Cummings prop- 
agators entering the phonon induced scattering terms. 
This leads to the interpretation that the phonons interact 
not with the bare electron, but rather with an electron- 
photon quasi-particle often referred to as a polariton. 
Indeed, if we approximate the U(t) matrix in the phonon 
induced scattering terms with the time-evolution opera- 
tor obtained for g = 0, i.e., the non- interacting QD-cavity 
system, then U(t) becomes strictly diagonal ■ . As a con- 
sequence g%(t) = and (i) = hr" 2 J* dtf D<(t') and 
the phonon induced scattering terms would not depend 
on the properties of the QD-cavity system. 



E. Polar on frame 



In the RDM formalism we derive an EOM for 
p(t)=Tr ph { X (t)}, 



(39) 



which is useful for calculating expectation values pro- 
vided that the operator of interest belongs to the system 
part of the Hilbert space. In this case we may perform 
the following operation 

(0(t)) = Tr s+ph { X (t)0} = Tr s {Tr ph { X (t)} 0} (40) 
= Tr s {p(t)0} . (41) 

If we now perform an arbitrary basis change operation 
given by the unitary operator T, where T^T = T~ 1 T — 
J, the expectation value of the operator O must of course 
not change, hence 

(0(t)) = Tr s+ph { X (t)0} (42) 

= Tr s+ph {TTt x (t)TTtOTTt} (43) 

= Tr s+ph {x(t)0} , (44) 

where the bar signifies the operator in the new basis. In 
the new basis we may also define a RDM for the system 
as follows 



p(i)=Tr ph {*(*)}. 



(45) 



However, in order for this object to be useful for calcu- 
lating physical expectation values, we need to be able to 
perform the following operation 

(0(t)> = Tr s+ph {x(t)0} = Tr s {Tr ph {*(*)} 0} (46) 
= Tr s {p(t)0} . (47) 



That is, the basis change should not entangle the system 
operator with the reservoir degrees of freedom or more 
formally O — o s <8> 7 p h, I p h being the identity operator in 
the phonon Hilbert space. 



In the case of the polaron transformation, see Eq. (Ill, 



all system projection operators are left invariant under 
the polaron transformation, i.e. a nm = cx nm , except for 
the off-diagonal operators: a\2, 013, and their hermi- 
tian conjugates. This has the consequence, e.g., that the 
bare electron polarization (c|(i)c g (t)) = Tr s [p(t)ai3], of- 
ten used to calculate the linear optical susceptibility, can- 
not be determined directly within the polaron framed 
Fortunately, all operators needed for our purposes are left 
invariant. 

As the polaron transformed Hamiltonian derived in 



Section II C is expressed in terms of bare QD-cavity op- 
erators, the elements of the RDM that are projected out 
are with respect to the bare QD-cavity system operators 
and hence do not always correspond to the actual physi- 
cal elements. To distinguish between expectation values 
calculated in the polaron and original frame, we intro- 
duce the following notation for the expectation values in 
the polaron frame 

(0(t)) p = Tv s {p(t)0}, (48) 

and as a consequence we get a new vector representation 
of the RDM in the polaron frame 



(aii(i))> 22 (i))> 12 (i)) D ,(a 21 (t)) t 



(tf 23 W)>3 2 W)>13W) p >3lW) p • (49) 



The polaron transformed Hamiltonian is given by 

H = H s i + -f/g'-ph' + #o,ph, (50) 



where the individual terms are defined in Eq. (14 1. As in 



the previous section we set up the EOM for the RDM 

d t p(t) = -ih~ x [H s ,,p(t)]+S LA (t) 

+ (L {a 32 ,K} + L {(731, T} + L {a ll7 2 7 }) p(t), (51) 

where the LA scattering term in this case contains the 
interaction Hamiltonian -ff S '_ p h'- The coupling matrices 
in the polaron frame for the coherent and Lindblad terms 
are identical to those in the original frame, see Eqs. ( p4"| , 
(25), and (26), except that the replacement g — > (X)g 



should be performed in the coherent terms. The terms 
arising from the coupling to the LA phonons are 

r-Tx(t) +r a (*) -iG* 2 (t) +iG 2 {ty 
+r 1 (t) -r 2 (t) +iG* 2 (t) -iG 2 {t) 

+iG 1 (t) +iGx(t) -7i(f) -iGt(t) 
_-iG\{t) -iG\{t) +iG 3 (t) -7?(i). 



(11) 



(t) 



and 



,(22) 



(*) = 



-72 (t) 



iG 4 (t) 






-7l(<) 


-iG\{t) 



iG 5 (t) 


-73 (t) 






-iGUt) 


"73 W. 



(52) 



(53) 



All elements are explicitly denned in Appendix [F] As 
these expressions are given in the polaron frame, we can 
not interpret the different terms as easily as in the orig- 
inal frame. However, we will still note a few differences 
and similarities. We now see a direct phonon induced 
lifetime renormalization of states 1 and 2 through Ti and 
r 2 , as well as several quantities playing a role similar to 
G^(t) in the original frame, via the G n (t)'s. Also, all 
polarizations now have a phonon induced pure dephas- 
ing rate, given by the quantities 7 n (t), associated with 
them. All quantities are composed from terms of the 
form 



K 
where 



nmki(t)=9 2 I dt'U n , m {t')UUt')B±{t'), (54) 



U(t) =exp(-iH a d/h), (55) 



is the time-evolution operator with respect to H s >. The 
functions B+(t ) are correlation functions for the polaron 
defined in Eq. (E15l and play a role similar to D<(t) in 
the original frame. The structure of K^ mU {t) is similar 
to that of the scattering terms in the original frame, but 
the interpretation is complicated by the fact that we are 
in the polaron frame. 



F. The long-time non-Markovian limit 

The scatterings terms arising from the TCL are time- 
dependent, giving rise to non-Markovian behavior. In the 
case of an initial excitation of the system, the duration of 
the time-dependence is set by the memory depth of the 
associated reservoir correlation function, D<(t) for the 
original and B± (t) for the polaron frame. This is evident 
from Eqs. pjj), p0|, p2l, and ((541) as the time-evolution 



operator itself for cither frame does not decay. 

In Fig. [2] we show examples of the various correla- 
tion functions for a range of relevant temperatures. The 
correlation function in the original frame, (i) , has a 



temperature independent imaginary part (see Eq. (|33j)), 
whereas the real part varies significantly with temper- 
ature. The amplitude is smallest and memory depth 
is largest for low temperatures (the memory depth is 
extracted from the normalized correlation function, not 
shown) , where an increasing temperature leads to a larger 
amplitude and smaller memory depth. In the polaron 
frame the corresponding correlation functions are B± (t) , 
for which both the real and imaginary part are temper- 
ature dependent. The amplitude and memory depth be- 
have as in the original frame. For completeness we also 



show ip(t) entering B±(t), see Eq. (E26) 



Above we discussed the dependence of the phonon cor- 
relation functions on temperature, however other param- 
eters also influence the amplitude and memory depth of 
the correlation functions. The spatial extent of the QD 
wavefunction turn out to be important. The phonon cou- 
pling matrix element, see Eq. ([6]), is directly related to 






FIG. 2. Illustrations of the various correlation functions 
for the phonon reservoir for fixed material parameters^ and 
a range of temperatures. Solid (dashed) lines are for the 
real (imaginary) part. D > (t) can be obtained as D > (t) = 
[D<(t))\ 



the spatial Fourier transform of the absolute square of 
the wavefunction of the relevant QD state. A small QD 
will have relatively wide spectrum in fc-space and thus 
couple to more phonon modes, causing the correspond- 
ing correlation function to decay faster. Conversely, a 
large QD will have a more narrow spectrum and couple 
to fewer phonon modes, resulting in a slower decay of the 
correlation function^. In the following, we keep the size 
of the QD fixed and will not investigate this further. 

From Fig. [2] we conclude that the time-dependence of 
the phonon correlation functions and therefore the TCL 
scattering terms only becomes important within the first 
few ps of the time evolution. For the time-dependence 
of the rates to have a significant effect on the dynamics, 
the RDM has to change significantly within the first few 
ps after the initial excitation, which is not the case for 
experimentally relevant parameters. For this reason, we 
may safely let t — > oo in all TCL scattering terms render- 
ing them as constants. While the long-time limit is well 
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justified for studying population decay dynamics, this is 
not the case for quantities depending sensitively on quan- 
tum coherence, e.g., the degree of indistinguishability of 
single photons^. 

Taking the t — > oo limit in the TCL is sometimes 
referred to as a Markov approximation^, whereas the 
non-Markovian regime is accessed for times smaller than 
the memory depth of the reservoir. In the case of a 
memory-less reservoir, the long time limit is exact and 
does not impose any further approximations. A memory- 
less reservoir is assumed in the derivation of the famous 
Lindblad result, see Eq. (18), which is customary re- 



ferred to as the Markovian limit in the field of cQED. 
In our model the reservoir does, however, have memory 
and we obtain qualitatively different results compared to 
a Markovian description of the phonon coupling within 
the Lindblad formalism, even though we take the long 
time limit in the TCL scattering terms. To distinguish 
the two qualitatively different descriptions, we will re- 
fer to the memory-less (Lindblad) case as the Markovian 
and the case including memory effects as non-Markovian, 
even though the t — > oo limit has been taken. 



IV. RESULTS 

In this section we present the results obtained from the 
theory described in the previous sections. In Section [lV A| 
we provide a parameter investigation of QD decay dy- 
namics obtained by numerically solving the EOMs in the 
time-domain and using the polaron frame. We chose the 
polaron frame in order to obtain the most accurate re- 
sults. In Section [lVB| wc derive analytical expressions for 
the QD decay rate within both the original and polaron 
frame. We compare them numerically and discuss the 
insights that are obtained from their analytical forms. 



A. Quantum dot decay dynamics 

In Fig. [3] we show a series of decay curves calculated 
within the polaron frame for an initially excited QD and 
compare the res ults for different signs and values of the 
detunin g^ 7 * 19 * 51 !. The excitation could be due to an op- 
tical pulse, resonant with the photon-emitting \g) O |e) 
transition or higher states of the QD. The chosen param- 
eter values (g > n, T, 7) places this system well within 
the so-called strong coupling regime and the temperature 
has been set to K to freeze out thermal excitation of 
phonons. 

For the resonant case we observe a very fast decay, 
and clear Rabi oscillations, indicating the strong cou- 
pling regime. For non-zero detuning we observe an 
asymmetry with respect to the sign of the detuning, 
which has be en pre dicted theoretically^ and observed 
experimentally^EQEU. The physical origin of the asym- 
metry is due to spontaneous emission of phonons, while 
absorption of phonons is unlikely at very low tempera- 




200 400 
Time, t [ps] 



600 



FIG. 3. QD decay curves for an initially excited QD, calcu- 
lated as <7ii(t). The curves are for different signs of the detun- 
ing, solid (dashed) is for negative (positive) detuning, denned 



as A = w e 



Parameters: T = K, hg — 150 peV, 



Hk = 100 ^eV, r = 1 ns _i , and h-y = ^eV. 



tures, which could otherwise restore symmetry. The de- 
cay is fastest for positive detuning, as here the initially 
excited electron may emit a phonon to become resonant 
with the cavity and decay through it, whereas for neg- 
ative detuning, the absorption of a phonon is required. 
It is clearly seen that the asymmetry is strongest for in- 
termediate detuning values, which may be explained by 
examining the interaction matrix element, sec Eq. ©. 
From the nature of the deformation potential interaction, 
the matrix element vanishes for small phonon energies 
becoming proportional to while for large energies 

the form factor imposed by the finite QD wavefunctiorP^l 
causes the matrix element to decay. This gives rise to a 
maximum in the phonon matrix element, leading to the 
largest degree of asymmetry. 

To more systematically quantify the dependence on 
detuning and the influence of finite temperature on the 
phonon induced asymmetry, we calculated the degree of 
asymmetry by taking the ratio between the slow QD life- 
time for A < 0, ta<o, and the faster lifetime obtained for 
A > 0, ta>o- The results are presented in Fig. [4] along 
with the absolute lifetime for both signs of the detuning. 
The lifetime is obtained by fitting a single exponential 
to the decay curve obtained from the numerical solution 
of the model. In the situations where the decay is os- 
cillatory the fitted lifetime thus represents the decaying 
envelope of the entire curve. 

For the zero temperature case studied in Fig.[3j we ob- 
serve a degree of asymmetry of almost 2 near a detuning 
of approximately 1 meV. The value of the detuning for 
which the maximum is obtained is determined by the ef- 
fective size of the QD through the form factor entering 
the phonon matrix element M$^. For comparison, we 
also show the curve with no phonons in the model and 
which shows that for low temperatures the QD lifetime 
for A < is only very weakly influenced by the phonons. 
As the temperature is increased the degree of asymme- 
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1 2 3 4 

Detuning, h\A\ [meV] 



FIG. 4. (Top) QD lifetimes for negative, ta<o, (solid curve) 
and positive, ta>o, (dashed curve) detuning at three temper- 
atures for a range detuning values. The black curve is with no 
phonons in the model. (Bottom) Degree of asymmetry quan- 
tified by the ratio between the QD lifetimes for opposite sign 
of detuning. Parameters are hg = 150 fieV, fire = 100 )J,eV, 
r = 1 ns~\ and fir/ = fj,eV. 
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FIG. 5. Dependence of (X), see Eq. (E14l, on temperature. 



try decreases. Intriguingly, the QD is seen to decay more 
slowly at very large detuning as temperature is increased, 
even though this is basically outside the bandwidth of the 
phonons. We believe this to be due to the renormaliza- 
tion of g caused by (X), lowering the effective value of g, 
see Fig. [5 where the temperature dependence of (X) is 
shown. The smaller asymmetry for higher temperatures 
is caused by the presence of thermally excited phonons, 
making it more probable for the electron to absorb a 
phonon and thereby becoming resonant with the cavity 
in the case when uj cav > w cg , i.e., A < 0. 

To illustrate the behavior of the phonons at different 
temperatures, we calculated the real part of the phonon 



correlation function Eq. ( 33 1 in the frequency domain 



Re [£>>(«;)] = 7r^|Af fe | 2 [n fc <5( W + u h ) 

k 

+ [n h + l]S(u-w h )], (56) 



where the Fourier transform is calculated as D > (lj) — 
L°° dte l ^ +M ^ t D > {t), where + is a positive infinitesi- 
mal. The quantity Ke[D > (uj)] gives information about 
the phonon modes interacting with the QD for a given 
temperature and can thus be considered as an effective 
phonon density. Also, it enters directly into the QD de- 
cay rate, as will be demonstrated in Section |iVB| 
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FIG. 6. Effective phonon density Re [D > (lj)], see Eq. (56 1, 
for a set of temperatures. The weak shoulder visible at low 
temperatures near hw ~ 2 meV arises due to different lo- 
calization lengths for the electron in the excited and ground 
states. 



In Fig. [6] we show Re [D > (co)] for a range of tempera- 
tures. For zero temperature, no phonons are available for 
absorption processes, corresponding to negative frequen- 
cies in the figure, while the vacuum phonon field reveals 
its presence through the non-zero density for positive en- 
ergies. This explains why the asymmetry is largest for 
zero temperature, as illustrated in Fig. [4] As the tem- 
perature is increased, more and more phonons are be- 
ing thermally excited and become available for both ab- 
sorption and stimulated emission processes. The strong 
asymmetry is no longer present in the effective phonon 
density, which correlates nicely with the observed behav- 
ior of the QD lifetimes. 

We will now investigate the dependence of the phonon- 
induced asymmetry on the light-matter coupling strength 
g. In Fig. [7] we show decay curves for a QD for both 
signs of the detuning and vary the light-matter coupling 
strength from very small values to large values represent- 
ing current state-of-the-art sampled^. The temperature 
is fixed at K. The first observation is the decrease of life- 
time for increasing g, consistent with the Purcell effecP^. 
Furthermore, we also observe an increasing asymmetry 
between lifetimes for positive and negative detuning val- 
ues as g is increased. This trend is seen more clearly 
in Fig. [8] where we show the degree of asymmetry as a 
function of detuning, for varying light-matter coupling 
strength g. It is apparent that one may go from a situa- 
tion of basically no asymmetry, obtained for a sample in 
the regime of weak or intermediate coupling strength^, 
to more than a factor of 2 in ratio between lifetimes in 
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FIG. 7. QD decay curves for an initially excited QD. The 
curves are for different values of the light-matter coupling 
strength, g, with the solid (dashed) lines being for a detuning 
of -1 (+1) meV. Parameters: T = K, M = 100 ^eV, V = 
1 ns _1 , and fir/ = /ieV. 
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Background decay rate, Y [ns _1 ] 

FIG. 9. Degree of asymmetry as a function of QD background 
decay rate. The detuning is fixed at ft|A| = 1 meV, while the 
light-matter coupling strength is varied. Other parameters 
are: T = K, fin = 100 /LteV, and fry = ^eV. 
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Detuning, h\A\ [meV] 

FIG. 8. As Fig. [4] except here the light-matter coupling 
strength is varied and T — K. 

state-of-the-art sampled^. This behavior might seem 
surprising at first, since, as independently of the value of 
the detuning, the electron has to emit a photon in or- 
der to decay to the ground state, regardless of whether a 
phonon was emitted or absorbed. From this observation 
one would expect the degree of asymmetry to be inde- 
pendent of g, since the Purcell enhancement scales with 
g, independently of the detuning. The reason for the 
dependence on g is simple, as will be explained below. 

The degree of asymmetry is seen to approach unity in 
the limit of small light-matter coupling strength, where 
cavity-mediated effects play a less significant role for the 
QD decay dynamics. Indeed, in the limit of small g or 
large A, the dominant decay channel for the QD becomes 
the background decay rate, T, which includes, e.g., decay 
into radiation modes and non-radiative decay. To illus- 
trate the effect of the background QD decay rate, we show 



in Fig. [9] the degree of asymmetry as a function of T for a 
few typical values of the light-matter coupling strength, 
covering weak, intermediate, and strong coupling. For a 
typical weak coupling sample, Kg — 30 /ieV, a noticeable 
asymmetry is only visible for very small V, corresponding 
to cavities where radiation modes are strongly suppressed 
such as photonic crystal cavities. The asymmetry disap- 
pears as the phonon contributions become dominated by 
the background decay rate. On the other hand, for a sam- 
ple well within the strong coupling regime, Kg = 150 ^eV, 
a significant asymmetry should be observable for basi- 
cally all values of the background decay rate. 



B. Approximate analytical expressions 

While the results from the previous section are nu- 
merically exact solutions for the dynamics, more physi- 
cal insight can be gained through approximate analytical 
expressions for the QD decay rates. In the limit of large 
detuning, A 3> g, such expressions can be obtained in 
both the original and polaron frame. This is possible as 
we can adiabatically eliminate the involved polarizations, 
and the time evolution operator, U (t) , may be expanded 
to a low order in the quantity g/A 7 see Appendix [G] for 
details. 

In the original frame we obtain the following expression 
for the total QD decay rate 

r to t = r 
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and for the polaron frame we obtain 



Ct = r 



> 7tot 

2 _j_ A 2 



7tot 



2g 2 Re[B-(uj = A)], (58) 



where the total dephasing rate is defined as 



1 



7tot = ^(k + T)+j. 



(59) 



In Eqs. ( |57| and ( 58 ) the Fourier transform is calculated 
as f(ui) = f™ dte i( - u+i0+ ^f(t), where 0+ is a positive 
infinitesimal. 
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FIG. 10. Comparison of QD lifetimes obtained through the 
appr oxim ate rates (dashed curves) in the original (blue), 
Eq. (57 1, and polaron frame (red), Eq. (581, and a single ex- 



ponential fit (solid curves) to the numerically exact solution. 
We also show the result when phonons are not included in the 
model (green). Parameters are: hn — 100 fieV, Kg — 150 /ieV, 
r = 1 ns , and tvy = ^ieV. 



In 



Figs. [10] and 11 we compare the QD lifetime 
(t = l/r to t) calculated from the approximate expres- 
sions with single exponential fits to the numerically ex- 
act solutions, for two typical sets of parameters. For all 
but very small detuning values, the approximate expres- 
sions compare very well to the corresponding numerical 
fits. The strong asymmetry at low temperatures, as well 
as the more symmetric decay rates at elevated tempera- 
tures, are well captured by the approximate expressions. 
At high temperatures, we observe significant deviation 
between the results in the original and the polaron frame. 
This is expected as only the polaron frame takes into ac- 
count multi-phonon effects that becom e increasingly im- 
portant at elevated temperature d 19 * 5 ^ . 

The expression for the decay rate in the original frame, 
Eq. (57), has a form very suitable for interpretation. In 
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FIG. 11. As in Fig. [10] except that Hk = 250 ^eV. 



addition to the background QD decay rate T, there are 
two contributions. The first contribution accounts for 
the direct decay of the QD through the cavity by emis- 
sion of a photon, with the total dephasing rate 7tot in- 
cluding a Lindblad pure dephasing rate t^h This gives 
rise to the familiar symmetric dependence on the detun- 
see the green curve in Figs. 



mg, see trie green curve m f igs. [10] and 11 However, 
the second contribution goes beyond the standard mod- 
els of cQED by depending on the effective phonon density 
Kc[D > (lj = A)] evaluated at the QD-cavity detuning, see 
Eq. (56) and Fig. M Thus, the phonon-assisted QD de- 
cay simultaneously depends on the cavity, through the 
Purcell rate prefactor, and on the availability of phonons 
that couple to the QD at the given QD-cavity detun- 
ing. Loosely, one can think of the second contribution 
as a product between the effective photon and phonon 
densities available for both spontaneous and stimulated 
processes. 

Based on the analytical expression for the QD decay 
rate in the original frame, Eq. ( 57 ) , we can provide a more 



physically transparent discussion of the dependence on g 
and r of the degree of asymmetry discussed in Figs. [7] 
[8] and [9] We begin by formally dividing the total QD 
decay rate into the three contributions discussed above 



Ttot 

r P = 2. 9 



r + r P + r ph , 

2 7tot 



7tot 



A 2: 



r 



P h 



7tot 



A 2 



Re [D>(u 



A) 



(60) 
(61) 

(62) 



where Y is the background decay rate, Tp is the usual 
Purcell enhanced rate, and r p h is the rate containing the 
phonon contribution and can be thought of as a phonon- 
assisted Purcell enhanced rate. With reference to Fig. [9] 
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FIG. 12. Contributions to the total QD decay rates given in 
Eq. (60 1. Note that for these figures a spherical QD model 



has been employed, using a harmonic confinement length of 
5 nm. Parameters: Hk = 100 /ieV, fry = peV, and T = K. 



we show in Fig. [12] the three contributions to r tot for 
two values, one small and one large, of T and g, as a 
function of detuning. For both values of the QD-cavity 
coupling, we observe that neither the bare Purcell rate 
nor the phonon-assisted rate are affected much by going 
from the small background decay rate, T — 0.2 ns -1 , 
to the larger background rate, T = 3.5 ns -1 . Close 
to resonance, also the total decay rate appears rather 
independent of the magnitude of the background as it 
is completely dominated by the bare Purcell enhanced 
rate. However, this picture changes dramatically once 
we increase the detuning and the contribution from the 
bare Purcell rate becomes comparable to the two other 
contributions. In the case of the large background rate 
and small QD-cavity coupling, kg = 30 [ieV, the con- 
stant background dominates over the phonon-assisted 
rate, Tph, and hardly any phonon-induced asymmetry 
is observed. Referring to Fig. [9] this situation corre- 



sponds to a typical micropillar cavity in the weak cou- 
pling regime. If we now decrease the background rate 
to a lower value, corresponding to a typical photonic 
crystal cavity in the weak coupling regime [Fig. [9j, 
the background and the phonon-assisted contributions 
become comparable and the degree of asymmetry con- 
sequently rises. This illustrates that one may enter a 
regime, where phonon-induced spectral asymmetries be- 
come significant, by changing the background decay, a 
parameter which is often thought of as being of minor 
importance and with trivial physical implications. In- 
creasing the QD-cavity coupling to values typically found 
in the strong coupling regime, hg — 150 /ieV, we signif- 
icantly increase both the bare and the phonon-assisted 
Purcell enhanced rates. For both values of the back- 
ground rate, a clear asymmetry in the total QD decay 
rate is now observed, owing to the fact that the con- 
stant and symmetric background rate no longer masks 
the phonon-assisted decay rates. 

The approximate expression in the polaron frame, see 
Eq. (58), is not as straightforward to interpret as the 
expression in the original frame. The background decay 
r enters in the same fashion and we also observe a term 
similar to the one representing decay directly through 
the cavity in the original frame. However, in contrast, 
the quantity (X) only enters the polaron frame, where it 
plays the role of renormalizing the light-matter coupling 
strength to a smaller value. The dependence of (X) on 
temperature is shown in Fig. [5] where it is seen that the 
renormalization can be quite significant. The last term 
involves the spectral properties of the phonons, through 
the Fourier transform of the correlation function _B_ (t) 



2g 2 Re[B-{uj = A)] = 
2g 2 (X) 2 Re 



dte iAt| eV (t) _A 



(63) 



where (p(t) is defined in Eq. (E22) and plays the role of a 
phonon-assisted QD decay rate analogous to Eq. (62) in 
the original frame. As B_(t) contains (X) as a factor, 
g is renormalized by (X) everywhere it appears. This is 
not the case for other cQED models also employing the 
polaron transformation™ The same formula has recently 
been independently derived and discussed by Roy and 
Hughes in Ref. [571 

The remaining part involving the Fourier integral over 
exp[ip(t)]~ 1 is harder to interpret than the corresponding 
expression for D > (uj) in the original frame. Even though 
ip(t) and D > (t) appear rather similar, compare Eq. (33) 
and Eq. (E22), Re[D > (u})] directly reflects the effective 
spectral features of the phonon reservoir. Also, in the 
original frame, D > (u) carries the familiar Lorentzian- 
style denominator of the cavity lineshape, which is miss- 
ing in the polaron frame. Mathematically, the Lorentzian 
denominator appears in the expression since the phonon 
induced term enters via a polarization, whereas in the po- 
laron frame, it enters directly as a lifetime. Despite the 
fact that they superficially look rather different, their 
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numerical values compare very well, especially for low 
temperatures, as evidenced in Figs. 10 and |ll| 



V. SUMMARY AND CONCLUSION 

In summary, we have presented a theory for cou- 
pled QD-cavity systems including the interaction with 
phonons and illustrated the importance of the phonon 
interaction for the QD decay dynamics. 

Furthermore, we have provid ed a d etailed account of 
the theory used in recent studies^- 11 * 20 !, which is based on 
a second order expansion in the phonon coupling, while 
accounting for the polaritonic nature of the QD-cavity 
to all orders. It was shown that it is essential to include 
the polaritonic nature in the interaction, when describing 
non-Markovian phonon reservoirs. 

For elevated temperatures, multi-phonon effects are ex- 
pected to play an important role. To study the influence 
of phonons in this regime, we included a theory based on 
the so-called polaron transformation, which takes certain 
phonon processes into account to infinite order, while still 
maintaining important polaritonic aspects of the QD- 
cavity system. 

Using the polaron theory, an extensive investigation 
of the parameter dependence of the QD decay dynamics 
was carried out for experimentally relevant regimes. An 
asymmetric detuning-dependence of the QD lifetime was 
observed, where a positive detuning, uj cg > w cav , yielded 
a significantly faster decay compared to negative detun- 
ing, uj cg < w cav . The faster decay observed for positive 
detuning reflects that the QD may emit a photon by the 
simultaneous emission of a phonon, thereby overcoming 
the energy mismatch. Conversely, for negative detun- 
ing, absorption of a phonon is required to bridge the gap 
in energy, but at low temperatures phonon absorption 
is very unlikely. As the temperature is increased, the 
asymmetry gradually disappears, due to the availability 
of phonon absorption processes. Apart from inducing 
spectral asymmetries, the interaction with phonons also 
gives rise to a significantly increased bandwidth of the 
QD-cavity interaction. It greatly extents the bandwidth 
beyond that imposed by the cavity linewidth normally 
thought to be the limiting factor, relaxing the resonant 
nature of many cQED phenomena. 

We also provide a simple explanation for the lack of 
experimental observations of phono n- induc ed asymme- 
tries in QD decay curves until recently ^ 11 * 51 !. We showed 
how the background decay rate of the QD, often consid- 
ered insignificant compared to other loss channels, plays 
a surprisingly important role in observing phonon effects 
for non-zero detuning. Phonon effects are strongest at 
relatively large detunings, 1 — 2 meV in our case, which 
typically spans many cavity linewidths of 0.05 — 0.3 meV, 
and thus the effect of the cavity is usually small at these 
detunings. In order for cavity-mediated effects, such as 
the phonon asymmetry, to remain significant either a 
small background decay or a large light-matter coupling 



strength is needed. Both of these requirements demand 
high quality samples, which have only become available 
recently. 

To provide further insight into the physics, we derived 
approximate analytical expressions for the total QD de- 
cay rate, which distills the essential ingredients added 
by the phonon interaction to well-known results from 
cQED. The power and accuracy of these expressions has 
recently been demonstrated experimentally and the effec- 
tive phonon density has been experimentally extracted^. 
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Appendix A: Equation of motion for reduced 
density matrix 

In this appendix we derive the equation of motion 
for the reduced density matrix of the QD-ca vity sy stem, 
which interacts with a large bosonic reservoi r 21 * 22 !. 

We start by defining the total Hamiltonian 



H{t) = H s (t) + H R + H SR = H (t) + H t 



Sli- 



(Al) 



where its(t) is the, possibly time-dependent, Hamilto- 
nian for the system of interest, H R is the Hamiltonian 
for the reservoir, and Hg R is the interaction between the 
two subsystems. For notational simplicity, we have in- 
troduced Ho(t) as the sum of the free contributions. 

The time evolution of the total density matrix, xWj 
is governed by the following equation in the Schrodinger 
picture 



ihd t x(t) = [H(t), x (t)] 



(A2) 



where H(t) is the Hamiltonian defined in Eq. (Al). We 
transform into the interaction picture with respect to 
Hs(t) + Hn, to facilitate a perturbation expansion in or- 
ders of the interaction fljsR. The transformation operator 
t/ff (t)(t, to) satisfies the Schrodinger equation 

ihd t U Ha(t) (t,t ) = {H s (t) + H R } U Ha(t) (t,to) 

= H (t)U Ho{t) (t,t ), (A3) 

where to is the initial time, and Uu (t){toi to) = /, with I 
being the identity operator. Ujj u-\(t, to) may be formally 
integrated, and due to the allowed time-dependence of 
the system Hamiltonian, we end up with the time-ordered 
expression 

£^o(t)(Mo)=T jexp^-ift- 1 ^ dt'H (t')^y (A4) 



15 



with T being the time-ordering operator. The interac- 
tion picture representation of the total density matrix is 
defined as 

X(t) = < (i) (Mo)x(^ff„ (i )(Mo), (A5) 
which leads to the following equation of motion for x(t) 



ihd t x(t) = [H sn (t),m 
This equation can be formally integrated 



X(t) = x(to) ~ ift 



dt' 



H S R(t'),x(t') 



(A6) 



(A7) 



By inserting this expression into the right hand side of 
Eq. ( A6 ) and tracing over the reservoir degrees of free- 



dom, we obtain a formally exact equation for the reduced 
density matrix of the system 



ihBtp(t) = Trn{[HsR{t),x(to)\} 



ilr 1 f dt'Tv K { \H SR (t), \H SR (t'), X(t')]]} , (A8) 

J to 



where p(t) = Tr R {x(*)} and Tr^ {• • •} denotes the trace 
operation with respect to the reservoir degrees of free- 
dom. To proceed further, we need to start invoking ap- 
proximations. The first approximation, known as the 
Born approximation, assumes that the total density ma- 
trix on the right hand side of Eq. ( A8 ) factorizes at all 



times, and especially at the initial time, hence 

X(t) = p(t)Ro & x(t) = p{t)R = p(t)R Q , (A9) 



where Rq — Rq is the density matrix for the reservoir, 
assumed to remain in a thermal state at all times and 
hence being time- independent. This approximation is ex- 
pected to hold for weak interaction between the system 
and reservoir. Motivated by the specific physical situa- 
tion considered, we shall assume that i?sR is written on 
the following form 



-ffsR = Pvv' B v 



(AlO) 



where P vl/ > is a pure system operator and B vv i is a pure 
reservoir operator. We assume that B vu > has the follow- 
ing property 



Tr R {R B VU ,} = (B vv ,) = 0. 



(All) 



If we now use eqs. ( A9 1 and (All ) in Eq. ( A8 ) we arrive 
at 



d t p(t) = -h~ 2 / di'Tr R { \H SR (t), \H SR (t'),p(t')Ro\\ } 

Jto 

(A12) 



which completes the formal derivation of the equation of 
motion for the reduced density matrix. 

To use the specific form of the interaction Hamiltonian, 
Eq. (AlO), we insert this into the above, expand the com- 



mutators, and rearrange the position of the B vv i , s with 
respect Rq to obtain well-defined expectation values over 
the reservoir operators. Performing these steps yields 



dtP(t) 



-fc" 2 / dt' £ { 



{p(t')P uW2 {t')P VlV2 {t) - P VlU2 {t)p(t')P v , v ,{t')} {B v[ „, 2 {t')B VlV2 {t)) }• (A13) 



In its present form Eq. ( A13 1 contains a memory in- 



tegral with p(t') as an integrand, therefore the time 
evolution depends on the past state of the system and 
therefore non-Markovian. However, it is well-known 
that a non-Markovian description may also be obtained 
in a fully time local theory, where the time evolution 
only depends on the present state of the system, but 
with time-dependent coefficients arising from the reser- 
voir interaction. One example of s uch a theor y is the 
timeconvolution-less approach (TCL^32 g | 35 | 53 | 58 l j n f ac t ; 
if one makes the replacement p(t') — > p(t) in Eq. (A13) 



the formal second order result in the TCL is obtained, 
which still describes a non-Markovian time-evolution. 
However, it is essential that this replacement is made 
within the interaction picture, where the only relevant 
time scale is the assumed slow time scale induced by the 
interaction with the reservoiiP^ We will present the re- 
sult for both the time-local and memory theory below. 



In Eq. ( A13 ) the time-evolution of the operators is only 



governed by the free Hamiltonian of the respective sub- 
system. Thus only the time-evolution operator for the 
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system 



U(t, t Q ) = T jexp (^-itr 1 jf dt'Hsit')^ | , (A14) 



should be used when transforming Eq. (A13) back to the 
Schrodinger picture. We obtain 



d tP (t) = -ih- 1 [H s (t),p(t)} + S(t), 



(A15) 



where we introduced the reservoir induced scattering 



term defined as 

S(t) = U(t > t )[d t p(t)]Ui(t > t ). (A16) 
By employing relations such as 

U(t,toW)U*(t,to) = U(t,t')p(l/)U\t,t'), (A17) 

and 

[/(Mo)iV 2 W f (Mo) = U(t,t')P VlV2 U\t,t'), (A18) 

we may derive the final form of the reservoir induced scat- 
tering S(t) term, defined in Eq. ( |A16 ), for both the time- 
local and memory version described above. The scatter- 
ing term with memory becomes 



•S'memW 



/ df E { 

•J tn II 

V\ V-2 v[ V 2 

{P vlV2 U{t,t')P v , iv!2P {t')U\t,t') - U(t,t')P u ^p(t')U\t,t')P Ul „ 2 } (B VlV2 (t)B <uk (t')) 

+ {U(t,t')p(t')P v , u ,tf(t,t')P VlV2 - P VlI/2 U(t,t')p(t')P v , v ,tf(t,t')} (B <v ,(t')B UxI/2 (t)) }, (A19) 



and the time-local one becomes 



StlW = -fr 2 J dt' £ { 

{P VlV2 U{tJ)P uW U\tJ)p{t) - U{t,U)P viv ,U\tJ) P {t)P VlV2 } (B„ lU2 (t)B u ,„,(t')) 

+ {p{t)U{t,t)P Kv 'U ] {tJ)P VlV2 - P UlU2 p(t)U(t,t')P v , v ,Ut(t,t')} (B v , v ,(t')B UlV2 (t)) }. (A20) 



Appendix B: The Hamiltonian 

In this appendix we describe the steps needed to ob- 
tain the Hamiltonian used in the main text, starting 
from a more fundamental Hamiltonian. The fundamental 
Hamiltonian is given by 

H = H s + i?o,ph + #c- P h + H 7 + H K + Hp- (Bl) 
The part governing the QD-cavity system is 

H s = huj g c^c g + fno c c\c c + htOc^a 1 a 

+ hg{^c\c c + clc s a). (B2) 
The free phonon Hamiltonian is 

-Ho.ph = ^2 huj k b{b k . (B3) 
k 



I 

The interaction between the electrons and the phonons 
is 

H e _ ph = {M^A c s + M eAc c ) (bl k + b k ). (B4) 

k 

The last three contributions iJ 7 , H K , and H-p refer to 
different reservoirs and their interaction with the system, 
giving rise to various forms of Markovian decay, which are 
introduced in the main text. Their explicit forms are not 
needed and will therefore not be discussed further in this 
appendix. For an elaboration on the above Hamiltonians, 
see the main text, Section |H A| 

We only consider a single electron in the system, hence 
the following relation holds 

c l c s + c t c c = 1, (B5) 
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which may be used to eliminate the ground state operator 
c|c g from the Hamiltonian. Physically, this elimination 
can be motivated by the fact that we only have one elec- 
tron in two levels, implying a perfect correlation between 
the two electronic states and hence it is sufficient to treat 
one of the levels explicitly. For reasons to be elaborated 
below, we choose the excited state. 

The elimination results in the following changes 

H s = hio cg clc Q + hu cav a' a + Hg(a i c g c e + CgC g a), (B6) 

where a constant energy term has been removed and we 
introduced the transition frequency of the QD, defined 
as w og = u c — ujg. Furthermore, we get 

ffc-ph = £ M k (bl k + b k )ctc c + ]T M* g (bl h + b k ), 

k 



k 



(B7) 



where we have introduced an effective interaction matrix 
element as M k = M k e — M k g and the last term without 
any electron operators corresponds to the phonon inter- 
action with the fully occupied ground state. In thermal 
equilibrium, before any excitation of the system, the QD 
is in its ground state and the phonon system is in an 
equilibrium state that takes into account the presence of 
the electron in the ground state. We wish to describe a 
situation that deviates from this thermal equilibrium and 
therefore it would be advantageous to take into account 
the phonon interaction with the fully occupied ground 
state from from very beginning^ This may be achieved by 
shifting the phonon operators^ through a unitary trans- 
formation, defined as 



b k -> e S b k e 5 , 



5' 



k fc 



(B8) 



This transformation only affects the phonon operators 
and leads to the substitution in the total Hamiltonian 



where B = £ fe M h (b f _ h + b k ) and a pq = \p) (q\. 

To simplify the equations, the QD-cavity detuning A 



can be introduced into Eq. (BIO). This can 



be achieved by moving into a rotating frame described 
by the unitary operator T(t) — exp(— icj cav [(7ii + <722]t), 
through which we obtain the Hamiltonian 



H s = hAau + fig(er 12 + cr 2 i). 



(B12) 



Appendix C: The polaron transformed Hamiltonian 

In this appendix we will apply the polaron transforma- 
tion to the total Hamiltonian. 

The total Hamiltonian presented in Appendix [B] con- 
tains contributions from reservoirs needed to include the 
Markovian losses. For the final equations, the decay 
terms arising from these Hamiltonians will, however, not 
be affected by the polaron transformation introduced in 
this appendix and they will be omitted in the follow- 
ing. We explicitly demonstrate this in Appendix [D] The 
Hamiltonian is 

H = ftA(7ii + %(<Ti2 + <T2l) 

+ a u ]T M k (bl k + b k )+J2 ^kb{b k . (CI) 

k k 

The transfor mation we apply is known as the polaron 
transformatiorP2EIH23] and is defined in the following way 



O = e"Oe 



where 



S = CTnC, 



C = J2 X Ubl k -b k ), \ k = F 



M k 



huj k ' 



(C2) 

(C3) 
(C4) 



For performing the transformation we employ the 
Baker-Campbell-Hausdorff formula which states 



b k -t b k 



M~ 



(B9) 



This removes the last term in Eq. ( B7 ) and introduces a 



new term given by — J2 k [2Af fc A/ gg K J / [huj k ] cjc c , which 
yields a simple energy renormalization that can be ab- 
sorbed into the bare excited state energy Huj e . 

For describing single photon emission, it is sufficient 
to operate in a one-excitation subspace of the QD-cavity 
Hilbert space. As a specific basis, we choose the follow- 
ing: {|1> = \e,n = 0) , |2) = |g,n = 1) , |3> = \g,n = 0)}. 
Along with the reformulations introduced above, project- 
ing the second quantized Hamiltonian on to this subspace 
changes the following parts of the total Hamiltonian de- 
fined in Eq. ETl 



H s = hu} cg an + hw cajV a 2 2 + %(ci2 + 021), (BIO) 

ffe-ph = -8*11, (Bll) 



O = e b Oe~ b 

= O + [S, O] + i [S, [S, O)] + i [5, [S, [S, O]]] 



(C5) 



The transformed operators are: 

' bk = b k - A_kcrii. (C6) 



*n = cm, 012 = 0i2e 



Inserting these expressions and simplifying the resulting 
Hamiltonian yields 

H = HAa n + hg(a 12 X + + a 2 \XJ) + £ hu k b\b k , 

k 

(C7) 

where the detuning has been redefined as 

A-J2\M k \ 2 /(h 2 w k ) (C8) 
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to take into account the so-called polaron shift of the |1) 
state and further we introduced the phonon operators 



(C9) 



While Eq. |C7| is still an exact representation of the 
original Hamiltonian, the electron-photon and electron- 
phonon interactions have now been mixed into a single 
term. One might say that the photons now interact with 
a polaron, the electron-phonon quasi-particle, instead of 
the bare electron. It would be advantageous to more 
clearly separate the electron-photon and the electron- 
phonon interaction. To achieve this separatio n 1 31 ! 32 ! we 
replace X± with X± + (X) - (X) in Eq. (|C7| to obtain 



H, 



0,ph, 



with 



H s > = HAa n + hg (X) (cr 12 + cr 2 i), 
H s >- P h> = hg(ai 2 SX + + cr 2 i<5A_), 

H , P h = y^^fe&fe&fe- 



(CIO) 

(Clla) 
(Cllb) 

(Cllc) 



where (X) is defined in Eq. (E14| and 5X± in Eq. (E16|. 
Now H s i contains what might be referred to as a system 
Hamiltonian, however, it is not the original system con- 
sisting of only the electron and photon, as the phonon 
quantity (X) has entered. It is, however, of great advan- 
tage to include this term in the new system Hamiltonian, 
since then photon processes are treated to all order as 
well as preserving the detailed balance conditional. This 
would not be case if the system Hamiltonian were defined 



as the first term in Eq. (C7), the reby ending up treating 
the photons only to second ordeaP^H The quantity (X) 
has the effect of renormalizing the light-matter coupling 
strength g. From its definition, Eq. (E14), it is clear that 



< (X) < 1, and hence the presence of the phonons will 
always decrease the effective light-matter coupling. The 
Hamiltonian H s >- r > contains the interaction between the 
system and reservoir, which has been made weaker by the 
introduction of the difference operators 5X± , making it 
more suitable for a treatment using perturbation theory. 



Appendix D: Lindblad decay rates under the 
polaron transformation 

In this appendix we will calculate the effect of the po- 
laron transformation on a typical Lindblad decay rate. 
We consider the radiative contribution to the background 
QD decay rate, which has complicated non-radiative con- 
tributions as well, which can not be treated in a simple 
manner. Our starting point is the Hamiltonian 

H = huj cg a cc + 2J hflia\ai + hu: k b\b k 

l k 

+ (7 o Mk ( b -k + b k )+Y^ %'(a|°-gc + <Weg), (Dl) 



describing a two-level QD with ground and excited states, 
{\g) , |e)}, coupled to a phonon bath given by the b k op- 
erators and a photon bath given by the a; operators. 
Applying the polaron transformation as described in Ap- 
pendix [C] we obtain 

H' = tUv' eg a ee + ^2 hQiajai + ^ huj k b k b k 

l k 

+ h 9i (a]x + <7 gc + ai X_cr cg ) (D2) 
i 

where uj' r(r includes the polaron shift and X± is defined in 



Eq. (C9). We now divide the transformed Hamiltonian 
as fol 



lows 

H' = #o + Hi, 
where the free part is 

H' = tiw' eg a ee + hVLia\ai + fku k b k b k , 

l k 

and the interaction part is 

H{ = '^2hgi(a\x + (T gc + a ; A_cr og ) 



= Ba 



(D3) 
(D4) 

(D5) 
(D6) 



where we have defined the combined photon-phonon op- 
erator B as 



B = J2hgiajx + = AX+. 



(D7) 



In the original frame the initial condition is assumed to 
be a fully factorized state 



X(0) = j OQD(0)®i?ph< 



photon ) 



(D8) 



where x(i) is the density matrix of the total system. Per- 
forming the polaron transformation on the initial density 
matrix entangles the QD and phonon operators, so that 
the initial state no longer remains fully factorized. This 
complicates the further application of the Reduced Den- 
sity Matrix formalism and is often neglected under the 
assumption that it is small 33 . Employing this approxi- 
mation we proceed with the following density matrix in 
the polaron frame 



X'(0)«PQD(0)(g>i?p h( 



-Rpho 



(D9) 



We now follow the standard procedure and can write 
down the EOM for the excited state population of the 



QD using Eq. (|A20 1 
d t n{t) = -tr 2 



d£ 



AH) 



(B(t-t')B^) +c.c. n(t) 
(D10) 
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From the assumption of a factorized density matrix we 
obtain 

e **«(t-f ) (bq. _ t ') B t) = (x+(t - t')X-) G(t - t'), 

(Dll) 

where the polaron correlation function (X + (t — t')X-) is 
given in Eq. ( E24[ ) and the photon correlation function is 



G(t-t') = Y J [H 



2 -i(n,- w ' )(t-t') 



(D12) 



If gi is approximately constant near Qi = uj' cg one has 

G(t-t') = h 2 TS(t-t'), (D13) 

where T is the photon-induced decay rate of the QD, 
while we neglect the photon Lamb shift. The equation 
for the QD decay now becomes 



d t n(t) 



-r / dt' 



S(t - t') {X + {t - t')X_) +c.c. n(t), 

(D14) 



where, due to the appearance of the delta function in the 
integrand we may use for the phonon correlation function 



(X+(t - t')X_) | i=i ,= (X+(0)X_> = 1, 



(D15) 



Therefor, within the stated approximations the polaron 
transformation docs not influence Lindblad decay rates. 



Appendix E: Properties of the phonon operators 

In this appendix we give various results related to the 
phonon operator arising from the polaron transformation 

Turk 

X ± =e ±c , C = Y> fe (&l fe -6 fe ), A fc = — . (El) 

k 

The operators X± ma y b e written in terms of so-called 
displacement operators^ 

D k (a) = exp (ab{ - a*6 fc ) . (E2) 

If we rewrite the operator C in the following way 

c = j2 x *( b -k - bk) = E( A -^i - A -^)< ( E3 ) 

k k 

we can write 

X ± = []cxp [±(\- k b{ - \*- k t>k)] = n^fc( ±A -fc)- 
k fe 

(E4) 

We will need the following useful propertie d 61 * 62 ! of the 
displacement operators 

D{{a) = D k 1 (a) = D k (-a), (E5) 
D k (a)D k (j3) = D k {a + 0) exp (flm [afi*]) , (E6) 
(Z? fc (a))=exp(-|a| 2 [n fc + l/2]). (E7) 



In the last expression 

(blb k ) 



n k 



exp (f3hu> k ) — 1 



(E8) 



is the average thermal occupation of phonons in mode 
k and j3 = (k^T)^ 1 is the inverse thermal energy. The 
brackets (• ■ ■) = Tr p h {p P h,o 1 1 '} denote the expectation 
value with respect to the thermal density operator for 
the phonons 



Pph,0 



exp(-/?g ph ,o) 
Tr ph {exp(-(3H phfi )} ' 



(E9) 



which can be written as a product of the density matrices 
for the individual k modes as p p h.o = life Pph c wnere 



Pphfi 



exp 



(-/3huj k blb k ) 



Tr ph , fc |exp(- J SSw fe 6|,6 fe )| 



(E10) 



is the density matrix for the feth phonon mode. 
The first property we will derive is 

(X ± (t)) = (X ± ) = (X) , 



(Ell) 



where the time-evolution is with respect to -ff p h,o> result- 
ing in the standard expression for free evolution 



bk(t) 



(E12) 



Combining Eqs. (E4) and ( |E12[) and t aking the thermal 
expectation value using eqs. |E9l and (E10) we getP^ 



(X±(t)) = [] W+e^A-fe)) , (E13) 

where the individual terms in the product may be evalu- 
ated using Eq. (|E7|, yielding 



(X±(t)) = []cxp (-|A fc | 2 [n k + 1/2]) 
fc 

= exp l-^|A fe | 2 [n fe + l/2] 



(X) 



(E14) 



being independent of time. 

Next we will evaluate the polaron correlation functions, 
defined as 



B+(t,t') = (8X±(t)5X±(t')) , (E15a) 
B_{t,t') = (6X± {t)SX T (<')) , (E15b) 

SX ± (t) =X±(t)- (X). (E16) 



where 



Inserting this into the definitions of B±(t — t') we easily 
find 



(5X a (t)8X b (t')) = ((X a (t) (X))(X b (t) (X))) 



(X a (t)X b (t')} - (X) 2 , 



(E17) 
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indicating that we only need to evaluate (X a (t)X b (t')) , 
where a, b = ±. From the above we get 



X a {t)X b (t') =J{D k {ae i ^ t \_ k )D k (be i ^ t ' \_ k ), 



and using Eq. (|E6j) allows us to write 

X a (t)X b (t') = n^fe(A- fe [ae iWfe< + be^ t '}) (E19) 

k 

x cxp [a6i|A fc | 2 sin(w fe (t - t'))] . (E20) 

Taking the thermal average and employing Eq. (E7) 
yields 

(X a (t)X b (t')) = cxp |-^|A fc | 2 (2n fc + 1) j 

x exp{-ab(p(t-t')}, (E21) 



(E18) where we have defined the funct 



ion 



ip(t - t') = ^|A fc | 2 {[2n k + 1] cos(tu k [t - t'}) -ism{uj k [t - t'])} 

k 

= E' Afc | 2 {^e l ^ [t ' t,] + [n k + i] e — I*-*'] } 



(E22) 
(E23) 



Comparing Eqs. ( |E21| ), ( |E14[ ), and ( |E22| we see that 



(X a (t)X b (t')) = (X) 2 e- ab ^ t ~ t '\ 



(X) 



(E24) 
(E25) 



Going back to Eq. (E15) and using Eq. (E24), we obtain 
the final result 



B ± (t,t') = B ± (t-t') = (X) 2 (e 



1 



(E26) 



where, as expected, the equilibrium phonon correlation 
functions depend only on the time-difference and not the 
absolute time. We will also be needing B±(t',t), i.e., 
with the time arguments interchanged. These functions 
are available through complex conjugation 

(5X a (t)6X b (t'))* = 
[Tr ph { Pphfi 6X a (t)5X b (t')}]* = 

Tr ph {[ PpK0 5X a (t)5X b {t')] f } = 

Tv ph {SX b {t')SX a (t) Pphfl } = 
(SX b (t')SX a (t)) , 

where the bar signifies multiplication by — 1 and we used 
(X + y = X-. Finally, we obtain the following relation 



B ± (t',t) = [B ± {t,t>)X 



(E27) 



Appendix F: Scattering rates in the polaron frame 

In this appendix we explicitly define the scattering 
rates in the polaron frame entering in Eq. ( 52 ) . The 



building blocks are the integrals K^ mkl (t) defined in 



. => ° nip v > 

Eq. (54 1. Viz the discussion in Section |III F we take the 
long-time limit in the integrals, hence we define 



ikl 



ikl 



(oo). 



The rates are 



r 2 

7i 

72 
73 

id 
iG 2 
tG 3 
iG 4 
iG 5 



2Re [K^ u + K+ 12 ] , 
2Re [K n22 + X+jj , 



K 



1122 

k: 



K, 



2211J ~ 
-^122l] 



K 



1221 



K, 



2112J 



L 1122j 
[^221l]*+[^2 + ll 



K. 
K 
K 



2122 



1112 



1221 



2221J 



2221J 



2221J 



2112J 

[-^2122 
[^1211 



K. 
K 
K 



2221 



1211 



1122 



2122J 



2122J 



2211 J 



(Fl) 



(F2) 
(F3) 
(F4) 
(F5) 
(F6) 
(F7) 
(F8) 
(F9) 
(F10) 
(Fll) 



Appendix G: Analytical expression for QD decay 
rates 



In this appendix we derive the analytical expressions 
for the QD decay rates discussed in Section [lVB We pro- 
ceed in two steps; Firstly, an expression is derived that 
is valid whenever the cavity can be adiabatically elim- 
inated, and secondly, we take the large detuning limit, 
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which simplifies the phonon induced rates to the expres- 
sion presented in the main text. We only perform the 
explicit derivation for the QD decay rate in the original 
frame, Eq. (57), but the derivation for the same quantity 
in the polaron frame, Eq. (58), follows a similar proce- 
dure. 

From |IIID| we get the EOM for the excited QD popu- 
lation 

dt (011(f)) = -r (011(f)) + 2«?Im [(012(f))] , (Gl) 
and the photon-assisted polarization 



where the total QD decay rate is defined as 



r + 2 5 - 



712 



7f 2 + A 2 



g + Re[G' 



712 



(G5) 



We are interested in the spontaneous emission rate from 
the QD due to the coupling to the cavity field, hence we 
expect the final result to scale with g 2 . For this reason 
we expand the QD-cavity evolution operator U(t) up to 
first order in g/A 



d t (a 12 (t)) 



-<A + 7i 9 ] (012(f)) 



-i[g + Q>] (011(f)) + 1 [g + G<] (022(f)) , (G2) 

with 712 = 7 + Re [712] + (k + T)/2 and where the long- 
time limit has been taken in all phonon-induced rates, 
wherefor we omit the time argument. For the cavity to 
be adiabatically eliminated it can not perform any back- 
action on the QD, hence it can not enter in the above 
EOM for the photon-assisted polarization and we put the 
cavity population, ((722(f)), equal to zero. This is valid in 
the regime where the cavity decay rate, k, is much larger 
than all other parameters. Furthermore, when the total 
dephasing time I/I712I is much shorter than the charac- 
teristic timescale for (011(f)), we may put d t (012(f)) = 0. 
From this we get 



(012(f)) = -i 



g + Q > 

—iA + 712 



which when inserted in Eq. (Gl ) yields 



d t (011(f)) = -r tot (^^f)) 



(G3) 



(G4) 



iAt, 



[/(f) = e """0-11 + 022 



1 » t - 

+ A (e 



l)(0i2 + 02i)+O((5/A) 2 ). (G6) 



Using this expansion and Eq. ( 29 ) we find 



Re[G>] oc 



Im[D>(u> = 0) - D>(uj = A)} 
A ' 



(G7) 



which is small compared to the remaining terms and will 
be neglected. From the expansion of the time evolution 
operator we also find that Re [712] scales as g 2 , which 
makes it a higher order effect that can be neglected. We 
finally arrive at Eq. (57) 



r tQt = r + 2. 9 2 2 lto ' 

7tot + A 2 



1 + — Re[D>(o; = A)] 
7tot 



(G8) 



where 7 to t = 7 + (k + T)/2 is the total dephasing rate. 
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